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Abstract 

We present the first evidence for the production of T(1D) states in the four-photon cascade, 
T(35) - 1X b(2P), Xb(2P) - 7 T(1D), T(LD) - -y X b(lP), Xb(lP) - 7^(15), followed by the 
T(15) annihilation into e + e~ or u + fi~. The signal has a significance of 10.2 standard deviations. 
The measured product branching ratio for these five decays, (2.5 ± 0.5 ± 0.5) • 10™ 5 , is consistent 
with the theoretical estimates. The data are dominated by the production of one Y(LD) state 
consistent with the J = 2 assignment. Its mass is determined to be (10161.1 ± 0.6 ± 1.6) MeV, 
which is consistent with the predictions from potential models and lattice QCD calculations. 

We also searched for T(3S) -» 7Xb(2P), Xb{2P) jT(lD), followed by either T(1D) -» ??T(1S') 
or T(1Z?) — > 7r + 7r~T(lS'). We find no evidence for such decays and set upper limits on the product 
branching ratios. 

PACS numbers: 14.40. Gx, 13.20.Gd 
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Long-lived bb states are especially well suited for testing lattice QCD calculations [l[ and 
effective theories of strong interactions, such as potential models |2j or NRQCD |3j. The 
narrow triplet-^ states, T(IS'), T(2S) and T(3S), were discovered in 1977 in proton-nucleus 
collisions at Fermilab [4|. Later, they were better resolved and studied at various e + e~ 
storage rings. Six triplet-P states, x&(2fj) and x&Q-Pj) with J = 2, 1, 0, were discovered in 
radiative decays of the T(3S) and T(2S) in 1982 p and 1983 jfj, respectively. There have 
been no observations of new narrow bb states since then, despite the large number of such 
states predicted below the open flavor threshold. 

In this paper, we present the first observation of the T(1Z?) states. They are produced in 
a two-photon cascade starting from the T(3S) resonance: T(35) — > 7Xb(2Pj), Xb(2Pj) — > 
7T(1D). To suppress photon backgrounds from 7r°s, which are copiously produced in gluonic 
annihilation of the bb states, we select events with two more subsequent photon transitions, 
T(1D) — > 7x&(l-Pj), Xfe(l-fj) — * 7Y(1<S'), followed by the T(15) annihilation into either 
e + e~ or fi + ^ (see Fig. [IJ. The product branching ratio for these five decays summing over 
T(lD 12i 3) contributions was predicted by Godfrey and Rosner |7J to be 3.76 ■ 1CT 5 . 

The data set consists of 5.8- 10 6 T(3S) decays observed with the CLEO III detector at the 
Cornell Electron Storage Ring (CESR). Charged particle tracking is done by a 47-layer drift 
chamber and a four-layer silicon tracker which reside in a 1.5T solenoidal magnetic field Q. 
Photons are detected using an electromagnetic calorimeter consisting of about 8000 CsI(Tl) 
crystals 0. The particle-identification capabilities of the CLEO III detector jic} are not 
used in the present analysis. 

We select events with exactly four photons and two oppositely charged leptons. The 
leptons must have momenta of at least 3.75 GeV. We distinguish between electrons and 
muons by their energy deposition in the calorimeter. Electrons must have a high ratio 
of energy observed in the calorimeter to the momentum measured in the tracking system 
(E/p > 0.7). Muons are identified as minimum ionizing particles, and required to leave 
150 — 550 MeV of energy in the calorimeter. Stricter muon identification does not reduce 
background in the final sample, since all significant background sources contain muons. Each 
photon must have at least 60 MeV of energy. We also ignore all photons below 180 MeV 
in the calorimeter region closest to the beam because of the spurious photons generated 
by beam-related backgrounds. The total momentum of all photons and leptons in each 
event must be balanced to within 300 MeV. The invariant mass of the two leptons must be 
consistent with the T(IS') mass within ±300 MeV. 

Much better identification of the T(IS') resonance is obtained by measuring the mass of 
the system recoiling against the four photons. The average resolution of the recoil mass is 
17 MeV. The measured recoil mass is required to be within —4 and +3 standard deviations 
from the T(IS') mass. The mass resolution of the produced T(1D) state depends on the 
measurement of the energies of the two lowest energy photons in the event. Thus, we require 
that at least one of them is detected in the barrel part of the calorimeter, where the energy 
resolution is best. The selected events are dominated at this point by T(3S) — > 7r°7r°T(15') 
transitions, which have a branching ratio an order of magnitude higher than the expected 
signal rate. In fact, the branching ratio measured for a subsample of events in which two 7r° 



candidates can be formed is consistent with the previous measurements [11|. To suppress 
this background, we require the invariant mass for any photon pair to be at least 2 standard 
deviations away from the nominal 7r° mass. The sum of such two tt° mass deviations squared 
must be larger than 6 for any pairing scheme. 

To look for Y(1.D) events, we constrain events to be consistent with a photon cascade 
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from the T(3S f ) to the T(IS') via one of the Xb(2Pj) and one of the Xb(lPj) states. Only 
J = 1 or 2 are used since the J = states have small decay fractions for electromagnetic 
transitions. For each J2P, J\p combination we calculate a chi-squared: 



2 / w \ \^ ( E ~ii ~ E "fi (M T (id), hp, Jip) \ 

XlD,J 2P ,J 1P { M T(lD)) = 2^ > 

3=1 V E ^ ) 

where E^j are the measured photon energies; £^ ec * e are the expected photon energies 
calculated from the known masses of the bb states and the measured photon directions in 
each event. The masses of the T(1D) states are not known. Therefore, we minimize the 
above chi-squared with respect to Mr (id) which is allowed to vary for each event. The 
above formalism requires that we know how to order the four photons in the cascade. While 
the highest energy photon must be due to the fourth transition, and the second highest 
energy photon must be due to the third transition, there is sometimes an ambiguity in the 
assignment of the two lower energy photons from the first two transitions, since the range 
of photon energies in the T(35) — > 7X6 (2-Pj) decay overlaps the similar energy range in 
the Xb{2Pj) ~^ 7^(1-0) transition. We choose the combination that minimizes the above 
chi-squared. There are four possible combinations of J2P, J\p values. We try all of them 
and choose the one that produces the smallest chi-squared, Xid = m ^ n Xid j 2P , j 1p - 

In addition to the four-photon cascade via the Y(LD) states, our data contain events 
with the four-photon cascade via the T(25) state: T(3S) — > 7Xfe(2Pj), x&(2Pj) — > jT(2S), 
T{2S) -> jXbO-Pj), XbilPj) -> 7T(15), T(15) -> l + l~ (see Fig.[TJ). The product branching 
ratio for this entire decay sequence (including T(1S) — > l + l~) is predicted by Godfrey and 
Rosner 0] to be 3.84 • 10~ 5 , thus comparable to the predicted Y(1D) production rate. In 
these events, the second highest energy photon is due to the second photon transition (see 
Fig. Q). Unfortunately, these events can sometimes be confused with the T(LD) events due 
to our limited experimental energy resolution. The second and third photon transitions in 
the T(25) cascade sequence can be mistaken for the third and second transitions in the 
1{1D) cascade sequence, respectively. Therefore, it is important to suppress the T(2S) 
cascades. We achieve this by finding the J2P, Jip (= 0, 1 or 2) combination that minimizes 
the associated chi-squared for the T(2S) hypothesis, x\s = m ^ n X2S j 2P Jip> wri ere xls ^ s 
exactly analogous to Xid with the Mr (id) replaced with Mr(2S)- We then require xls > 12. 
Notice that the masses of all intermediate states are known for the T(2S*) cascade, thus this 
variable is more constraining than Xid- 

To further suppress the T(2S) cascade events, we construct a quasi- chi-squared variable, 
X25, that sums in quadrature only positive deviations of the measured photon energies 
from their expected values. This variable is less sensitive than x\s ^° fluctuations in the 
longitudinal and transverse energy leakage in photon showers that sometimes produce large 
negative energy deviations and correspondingly a large x\s value. With the additional 
criteria xls > 3 and Xid < 10> the cross-feed efficiency for T(2S') events is reduced to 
0.3%, while the signal efficiency is 12%. The 7r°7r° background cross-feed efficiency is 0.02%. 
Monte Carlo simulation of the signal events is based on the photon transition rate predicted 
for the .1 = 2 T(1D) state by Godfrey and Rosner 0. We use the J = 1 assumption to 
estimate the model dependence of the signal efficiency. The proper angular distribution of 
the first photon in the cascade, T(3S') — > 7Xb(2P), is taken into account, resulting in a 4% 
relative change of the efficiency compared to the uniform distribution. Angular correlations 
in the subsequent photon transitions are neglected. 
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The data x\d distribution after all these cuts is shown by the solid histogram in Fig. EH- 
A narrow peak near zero is observed, just as expected for T(1D) events. The signal Monte 
Carlo distribution for T(1D) events is shown by the solid histogram in Fig. Eb- The back- 
ground Monte Carlo distribution for the T(2S) cascades, after a factor of 10 enhancement 
relative to the T(1D) normalization, is also shown for comparison. The T(3S) — > 7r 7r°T(15') 
Monte Carlo distribution is shown without the ir° veto cuts to increase the statistics. We 
conclude that the backgrounds cannot produce as narrow a peak as observed in the data. 

After all the selection cuts, we observe 38 events in the data with Xid < 10- The 
background estimates are 1.5 ± 1.4 and 1.3 ± 0.9 T(25) and 7r°7r° events, respectively. The 
errors on the background estimates include systematic effects. Feed-across from the other 
photon and hadronic transitions is found to be negligible. Continuum backgrounds, for 
example due to radiative Bhabha scattering events, were estimated to contribute 0.7 ± 0.7 
events, using data taken at the T(IS') resonance. After the background subtractions, the 
estimated signal yield is 34.5 ± 6.4 events. 

An alternative background subtraction method is obtained by fitting the Xid distribu- 
tion in the range between and 100 to the Monte Carlo predicted signal and background 
contributions. In this method the background normalization is effectively determined by the 
event yield observed in the tail of the x\d distribution. The background shape is assumed 
to follow the 7r°7r° Monte Carlo distribution with the tt° veto cuts removed to increase the 
Monte Carlo statistics (see Fig. |2t>) - A linear background fit was also tried and yielded 
similar results. The T(2S) background is fixed in this fit to the Monte Carlo simulation, 
normalized to the rate predicted by Godfrey and Rosner, since unlike all other backgrounds 
it tends to peak near the signal region. This method yields 38.5 ± 6.8 signal events with a 
signal efficiency of 13% in the extended Xid range. 

The significance of the signal is evaluated from the change of likelihood between the 
nominal fit and when fitting the data with the background shapes alone and corresponds 
to 10.2 standard deviations (8.9cr for 7777/i + /i~ and 5.1cr for 7777e + e~). The signal prod- 
uct branching ratio obtained with both methods of background subtraction is the same, 
£>(7777Z + Z~)y(id) = (2.5 ± 0.5 ± 0.5) ■ 10 -5 . Throughout this paper we quote branching 
ratios averaged over the and e + e~ channels. The first error is statistical, while the 

second error is systematic. The systematic error includes uncertainty in the background 
subtraction (8%), model dependence of the efficiency (8%), uncertainty in the detector sim- 
ulation (8%) and the number of T(3S) decays (2%). This branching ratio is consistent with 
the theoretically estimated rate 0]. 

A straightforward way to measure the mass of the produced T(1D) state is to calculate 
the mass of the system recoiling against the two lower energy photons in the event. This 
distribution is shown in Fig. EH- The width of the observed peak is consistent with the 
detector resolution, implying the data are dominated by production of just one T(1D) 
state. We use the signal line shape obtained from the Monte Carlo simulations to fit the 
data and determine the mass of this state. 

Another estimate of the true T(1D) mass is given by the mass value that minimizes 
Xid- This distribution is shown in Fig. Ek>- The data are again consistent with the single- 
peak hypothesis. The fit to the expected signal shape from Monte Carlo simulations is 
superimposed in the figure. While this method has a mass resolution of about 3 MeV, 
compared to a value of about 7 MeV for the missing-mass technique, the signal shape here 
has a complicated tail structure originating from photon energy fluctuations which can make 
a wrong J 2 p, J\p combination produce the smallest chi-squared value. This produces small 
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satellite peaks on both sides of the main peak. Both methods of mass determination give 
consistent results. The mass of the observed state is determined to be (10161.1 ± 0.6 ± 1.6) 
MeV, where the first error is statistical and the second systematic. The systematic error 
includes the measurement method dependence (±1.2 MeV) and the mass calibration error 
(±1.1 MeV). The significance of a possible second peak around 10175 MeV is only 1.9 
standard deviations. The recoil mass distribution discussed in the previous paragraph and 
shown in Fig. |Hk favors no events in any second peak. Thus, we cannot claim observation 
of more than one T(1D) state based on our data. 

The Xid minimization favors the J 2 p = 1, J\p = 1 cascade path for most of the ob- 
served events, indicating that the observed state is either J\r> = 1 or 2. Theoretically, the 
production rate of the J\d = 2 state is expected to be 6 times larger than for the J±d = 1 
state |7|. Therefore, we interpret our signal as coming predominantly from the production 
of the T(1 3 D 2 ). Small contributions of the T(l 3 Di) and Y(l 3 -D 3 ) with masses close to the 
observed T(1 3 D 2 ) mass cannot be ruled out. However, they are impossible to quantify from 
our data alone without prior knowledge of the fine-structure mass splitting. 

The measured mass is in good agreement with the mass of the T(1 3 D 2 ) state predicted 
by lattice QCD calculations |l| and those potential models which also give a good fit to the 
other known bb states 12] • All potential model calculations predict the T(1 3 D2) mass to 



be between 0.5 and 1.0 MeV lower than the center-of-gravity (c.o.g.) mass for this triplet. 
Adding this theoretical input to our results, we obtain (10162 ± 2) MeV for the c.o.g. mass, 
where we assigned an additional uncertainty of 1 MeV to the correction for the 1 3 -D2 — c.o.g. 
mass difference. 

Voloshin recently suggested that the r\ transition could be enhanced in T(1D) — > Y(IS') 
decays k|. Since the rj often decays to two photons, we can look for it in the same sample 
preselected for the four-photon cascade analysis. We reverse the Xid cu ^ (Xid > 10) to 
suppress the four-photon cascades via the T(1D) states. Otherwise they would contribute 
a smooth background to our rj search variable (defined below). Since we still want the two- 
photon cascade to produce a D— state via T(3S) — > 7x&(2P 2 ,:l), Xfc(2P2,i) — ► 7Y(1.D), we 
require that one of the two lowest energy photons fits the T(35) — > 7Xf,(2P 2 ,i) transition 
(70.0 < P 7 < 110 MeV). Because the backgrounds are small, we did not constrain the 
second photon energy and therefore we did not restrict the sample to any particular value 
of T(1P) mass. The signal efficiency is 13% (not including B(r] — > 77)). To search for 
the eta we analyze the invariant mass distribution for the two most energetic photons. The 
distribution of (M 77 — M v )/<tm for the data is shown in Fig. HI where <tm is the expected rj 
mass resolution. No signal is observed. To estimate the upper limit we fit this distribution 
with the eta line shape and a smooth approximation for the background obtained from 
the Monte Carlo simulations. The corresponding 90% C.L. upper limit on the product 
branching ratio is: B(T(3S) -» jjT(1D))B(T(1D) -» T]T(1S))B(T(1S) -» < 0.6-10~ 5 
or B{T{3S) -> ^T{1D))B(T(1D) -> rjT(lS)) < 2.3 • 10" 4 if we use the world average 
value for B(T(1S) — > l + l~) jll|. A systematic error of 8.3% is included by scaling up the 
upper limit by one unit of the systematic error. Dividing the estimated upper limit by 
the measured product branching ratio for the four-photon cascade, we obtain: B(Y(1D) — > 
T]T(1S))/B(T(1D) — > 77T(15)) < 0.25 (at 90% C.L.). Common systematic errors were 
taken out in this calculation. 

Predictions for the branching ratio of T(1D) — > tt + tt~T(1S) vary by orders of magnitude 
among various theoretical predictions (from 0.2% to 49%) |l4|. To look for these transitions, 
we selected 777r + 7r~/ + /~ events using similar selection cuts to our 7777/+/" analysis. After 
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requiring the di-lepton mass and the recoil mass against the 777r + 7r~ to be consistent with 
the T(1S) mass, and checking that the total momentum of the event is consistent with zero, 
we require at least one photon to have an energy in the 70 — 110 MeV range, corresponding 
to the T(3S) — > 7x fe (2P 2 ,i) transition. We then measure the mass of the intermediate bb 
state, assuming that it is produced by the two-photon cascade. This mass can be estimated 
by using either the photons or the pions. To get the best estimate, we average the two mass 
estimates by giving them weights inversely proportional to the mass resolution squared, 
as determined by Monte Carlo simulations. The weights are 40% for the 77 recoil mass, 
and 60% for the mass obtained using 7r + 7r~. The signal efficiency is 19%. The resulting 
mass distribution is shown in Fig. The prominent peak observed in the data is due to 
T(3S) -> 7Xfe(2P), Xb(2P) -> 7T(25), T(2S) -»■ 7i + n~T(lS). From a fit to this peak, we 
determine the product branching ratio for this T(2S) decay signal to be 1.13 ± 0.16 times 
the value derived from the individually measured transition rates . This provides a good 
check for our detection efficiency. 

There is no indication of any excess of events at the T(1D) mass value observed 
in our four-photon cascade analysis. To estimate an upper limit on the signal rate, 
we fit the data with a signal fixed at our observed Y(1 3 D 2 ) mass and a smooth back- 
ground parameterized by a cubic polynomial. The following limits (90% C.L.) are ob- 
tained: B(T(35) -> 77T(LD 2 ))i3(r(l J D 2 ) -> n + n-T(lS))B(T(lS) -> l+l~) < 2.7 ■ lO" 6 
or B(T(35) -> 77T(1 J D 2 ))£(T(1 J D 2 ) -> 7r+7r-T(l£)) < 1.1 • 10~ 4 . Dividing our up- 
per limit by the measured rate for the four-photon cascade we obtain: B(T(1D 2 ) — > 
tt + tt-T{IS))/B(T{ID 2 ) -> 77T(15)) < 1.2 (at 90% C.L.). We also set an upper limit 
for the production of any T(1D) state (followed by ir + n~T(lS) decay) with a mass in the 
10140 — 10180 MeV range, which comfortably covers the predicted size of fine-structure split- 
ting for the T (ID) triplet 0- Here, we do not try to subtract backgrounds and accept all 9 
events observed in this mass range as signal candidates. This results in the following upper 
limits: B(T(35) -► ^T(1D))B{T(1Dj) 7r+7r-T(lS))B(T(LS) -» l + l~) < 6.6 • 10~ 6 or 
B(T{3S) -> ^T(1D))B{T{1Dj) -> n+7i-T{lS)) < 2.7 ■ 10" 4 for a sum over all different 
J id values. 

These upper limits are inconsistent (lower by a factor of about 7) with the rate estimated 
by Rosner 0] using the Kuang-Yan model for Y{T{1D) — > 7r + 7r~Y(LS3) ^| and a factor 
of about 3 higher than the predicted rate based on the model by Ko JlfiJ. Our upper limits 
are about 30 times higher than those predicted by Moxhay's model |l7j |. 

In summary, we present the first significant evidence for the production of the T(1D) 
states in the four-photon cascade T(3S) — > x&(2P) ~~ > T(1D) — > Xb(lP) —> T(1S). The 
data are dominated by the production of one T(1D) state, consistent with the J = 2 
assignment. Its mass is determined to be (10161.1 ± 0.6 ± 1.6) MeV, in agreement with 
the potential models and lattice QCD calculations. The measured product branching ratio, 
(2.5 ±0.5 ±0.5) -10 -5 , is consistent with the theoretical estimate, especially when comparing 
with the predicted rate for the T(1D 2 ) state alone, 2.6 • 10~ 5 [jj. 

We have also searched for T(3S f ) — > 7X&(2P), Xb(2P) 7^(1-D) followed by either 
T(1D) — > rjT(lS) or T(1D) — > 7r + 7r _ T(lS'). We find no evidence for such decays and set 
upper limits on the product branching ratios. The latter are inconsistent with the Kuang- 
Yan model which predicts a large Y(LD) — > 7r + 7r~T(15') width. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, and the Texas Advanced 
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FIG. 1: The expected bb mass levels. The four-photon transition sequence from the T(3S') to the 
T(1S) via the Y(1D) states is shown (solid lines). An alternative route for the four-photon cascade 
via the T(2S) state is also displayed (dashed lines). 
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FIG. 2: Distributions of xId f° r ( a ) data and (b) Monte Carlo simulations of the signal and 
backgrounds. The solid histogram in (a) represents the data, while the dashed line represents the 
background fit described in the text. The solid histogram in (b) represents the T(1Z?2) signal 
Monte Carlo. The dashed histogram shows the simulated background from the T(2S) cascades. 
This distribution is scaled up by a factor of 10 in efficiency normalization to make it visible when 
superimposed on that of the signal Monte Carlo. The dotted histogram shows the Monte Carlo 
distribution for 7r°7r° transitions with the 7r° cuts removed, normalized to the number of entries in 
the T(2S) cascade background histogram. The vertical line indicates the cut value used for the 
Y(1D) mass analysis. 
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FIG. 3: Distributions of the measured Y(IZ)) mass in the data using (a) the recoil mass method, 
and (b) the x\d ^ method. The results of fits for a single T(1D) state are superimposed. The 
Xid fit method produces satellite peaks as explained in the text. 
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FIG. 4: Distribution of the deviation of the two-photon mass from the r/ mass divided by the 
estimated mass resolution for T(1D) — ► rjT(lS) candidates from the data (solid histogram) and 
from the signal Monte Carlo simulation (dashed histogram). 
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FIG. 5: The invariant mass distribution for the system recoiling against the two photons in 
T(3S) — > 777r + 7r~ Y(IS') events. The observed peak is due to transitions via the T(2S) state, 
followed by T(2S) — > 7r + 7r _ T(lS'). The arrow indicates where the signal due to transitions via the 
T(1 3 Z?2) state is expected. 
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